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The large-scale energy spectrum in tvifo-dimensional turbulence governed by the surface 
quasi-geostrophic (SQG) equation 

is studied. The nonlinear transfer of this system conserves the two quadratic quantities 
4-1 = ([(-A)i/4V]^)/2 and = {[i-A)^^^-)jj]'^)/2 (kinetic energy), where (•) denotes 
a spatial average. The energy density ^2 is bounded and its spectrum 5'2(fc) is shal- 
lower than in the inverse-transfer range. For bounded turbulence, 5*2 (fc) in the 
low-wavenumber region can be bounded by CA: where C is a constant independent of 
k but dependent on the domain size. Results from numerical simulations confirming the 
theoretical predictions are presented. 



1. Introduction 

The dynamics of a three-dimensional stratified rapidly rotating fluid is characterized by 
the geostrophic balance between the Coriolis force and pressure gradient. The nonlinear 
dynamics governed by the first-order departure from this linear balance is known as quasi- 
geostrophic dynamics and is inherently three-dimensional. The theory of quasi-geostrophy 
is interesting and the research performed on this subject constitutes a rich literature (see, 
for example, Charney 1948, 1971; Rhines 1979; Pedlosky 1987). This theory renders a 
variety of two-dimensional models that are appealing for their relative simplicity and yet 
sufficiently sophisticated to capture the underlying dynamics of geophysical fluids. One 
such model, the so-called surface quasi-geostrophic (SQG) equation, is the subject of the 
present study. 

Quasi-geostrophic flows can be described in terms of the geostrophic streamfunction 
ijj{x,t). The vertical dimension z is usually taken to be semi- infinite and the horizontal 
extent may be either bounded or unbounded. Normally, decay conditions are imposed 
as 2; ^ cxD. At the flat surface boundary z = 0, the vertical gradient of ■(/'(a;,t) matches 
the temperature fleld T{x,t), i.e. T{x,t)\z=o = dzip{x,t)\z=o. For flows with zero po- 
tential vorticity, this surface temperature field can be identified with {—Ay^'^^p, where 
A is the (horizontal) two-dimensional Laplacian. Here, the operator (— A)^/^ is defined 
by (— A)^/^'0(fc) = fc-0(fc), where k = \k\ is the wavenumber and V'(fc) is the Fourier 
transform of ■0(a;). The conservation equation governing the advection of the tempera- 
ture (— A)^/^-0 by the surface flow is (Blumen 1978; Pedlosky 1987; Pierrehumbert, Held 
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& Swanson 1994; Held et al. 1995) 

9t(-A)i/V + ^(^,(-A)i/V) = 0, (1.1) 

where J{ip, (f) = dx^dycf) — dyLpdx4>- This equation is known as the SQG equation. 

In this paper a forced-dissipative version of Hl.l() is studied. A dissipative term of the 
form fiAip, where /i > 0, which resuhs from Ekman pumping at the surface, is considered 
(Constantin 2002; Tran 2004). Since {—AY^^ip is the advected quantity, this physical 
dissipation mechanism corresponds to the (hypoviscous) dissipation operator A)^/^. 
The dissipation coefficient fi has the dimension of velocity and is not vanishingly small 
in the atmospheric context (Constantin 2002). The system is assumed to be driven by 
a forcing /, for which the spectral support is confined to wavenumbers fc > s > (in 
bounded turbulence, wavenumber zero is replaced by the minimum wavenumber). Thus, 
the forced-dissipative SQG equation can be written as 

at(-A)i/V + J(^, (-A)i/2V') = fiA^ + /. (1.2) 

It is customary in the classical theory of turbulence to consider a doubly periodic domain 
of size L; the unbounded case is obtained via the limit L — > oo. 
The Jacobian operator J(-, •) admits the identities 

(xJ(^,</')) = -(Mx,0)) = -(M^,x)), (1.3) 

where (•) denotes the spatial average. As a consequence, the nonlinear term in (|1.2I) obeys 
the conservation laws 

(VJ(V^, (-A)^/^)) - ((-A)i/2V.J(V', (-A)^/^)) = 0. (1.4) 

It follows that the two quadratic quantities ^'g = (|(— A)^/^'0p)/2 = / '^g{k)dk, where 
9 = 1,2, are conserved by nonlinear transfer. Here, 'i'e{k) is defined by '^e{k) = k^'^{k), 
^{k) is the power density of V' associated with wavenumber k and is a real number. 
Note that ^'2(fc) is the kinetic energy spectrum and ^'2 is the kinetic energy density. 

The simultaneous conservation of two quadratic quantities by advective nonlinearities 
is a common feature in incompressible fluid systems in two dimensions. Some familiar 
systems in this category are the Charney-Hasegawa-Mima equation (Hasegawa & Mima 
1978; Hasegawa, Maclennan & Kodama 1979) and the class of a turbulence equations 
(Pierrehumbert et al. 1994), which includes both the Navier-Stokes and the SQG equa- 
tions. These conservation laws, together with the scale-selectivity of the dissipation and 
unboundedness of the domain, are the building block of the classical dual-cascade theory 
(Fj0rtoft 1953; Kraichnan 1967, 1971; Leith 1968; Batchelor 1969). This theory, when ap- 
plied to the present case, implies that ^'i cascades to low wavenumbers (inverse cascade) 
and ^'2 cascades to high wavenumbers (direct cascade). For some recent discussion on the 
possibility of a dual cascade in various two-dimensional systems, including the Navier- 



Stokes and SQG equations, see Tran fc Shepherd (2002) Tran & Bowman (2003b,2004) 



and Tran (2004) The inverse cascade toward wavenumber A; = would eventually evade 
viscous dissipation altogether because the spectral dissipation rate vanishes as fc 0. 
Hence, according to the classical picture, necessarily grows unbounded, by a steady 
growth rate d^'i/di > 0, as t ^ 00. Strictly speaking, one may have to address the 
possibility of a dissipated inverse cascade, i.e. one for which the dissipation of ^'i occurs 
at scales much larger than the forcing scale and for which d^'i/dt has a zero time mean. 
Such a cascade is not a plausible scenario (and is not the traditional undissipated inverse 
cascade) in fluid systems, dissipated by a single viscous operator, where the viscous dis- 
sipation rate diminishes toward the large scales. A discussion of this issue can be found 



in Tran (2004) 
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In this study, upper bounds are derived for the time averages of the kinetic energy 
density \l/2 and of the large-scale spectrum 5'2(fc). These bounds are derived from the 
governing equation, involving simple but rigorous estimates. The bound on \E'2 is valid 
in both unbounded and bounded cases, and a straightforward consequence of this bound 
is a bound on the energy spectrum, which also applies to both unbounded and bounded 
turbulence. Another bound on the large-scale energy spectrum is derived by estimating 
the nonlinear triple-product term representing the inverse transfer of ^'i. This result ap- 
plies to bounded turbulence since upper bounds for the triple-product term are inherently 
domain-size dependent. The difficulties of extending this result to the unbounded case are 
discussed. Some numerical results confirming the theoretical predictions are presented. 



2. Bounded dynamical quantities 

A notable feature of unbounded incompressible fluid turbulence in two dimensions 
is the appearance of infinite quadratic quantities (per unit area): namely, the kinetic 
energy density for Navier-Stokes turbulence and '^i for the SQG case. According to 
the classical theory (applied to the SQG case), a (steady) injection of ^'i, applied around 
some finite wavenumber s, cascades to ever-larger scales, leading to an unbounded growth 
of (this is presumably the case for the general quadratic invariant V^q in the so-called 
a turbulence; c/. Tran 2004). In other words, if the classical inverse cascade is realizable, 
unbounded incompressible fluid turbulence in two dimensions constitutes an ill-posed 
problem, in the sense that a key quadratic invariant becomes infinite. Of course, there 
still exist finite quadratic quantities, in particular the dissipation agent for each quadratic 
invariant. This section is concerned with these quantities. 

On multiplying H1.2|l by and (— A)^/^^ and taking the spatial average of the resulting 
equations, noting from the conservation laws that the nonlinear terms identically vanish, 
one obtains evolution equations for ^E*! and ^'2, 

^^^, = -2n-^2 + {M). (2.1) 

Avl;^^ -2^**3 + (/(-A)i/V). (2.2) 
at 

Using the Cauchy-Schwarz and geometric-arithmetic mean inequalities, one obtains up- 
per bounds on the injection terms in 12.1|l and H2.2|l : 

(/(-A)i/V^)<(|(-A)3/VP)^/^(|(-A)-VVP)i/2</iVl/3 + ^-iF_i, (2.3) 

where the 'integration by parts' rule ((-A)'^(/)x) = {{- ^)''" x) , for = 9' + 6", 
has been used and Fe = (|(-A)'^/'*/|2)/2. Substituting in (jTIJ) and yields 

< -At*2 + M"'-F^-2, (2.4) 

at 

^*2 < -A^*3 + M"'i^-i. (2.5) 

To avoid unnecessary complications, zero initial conditions are assumed, so that for T > 
the time means (d\E'i/dt)t = 4'i(T)/T and (d^2/dt)t are non-negative. One can then 
deduce upper bounds on the time means (^2)* and (^3)f, which are valid regardless of 
whether or not remains finite in the limit t 00: 

{^2)t<ti-^{F^2)u (2.6) 
{'^i)t<l^-^{F-i)t. (2.7) 
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For 6 G (2,3), (^'e)* is also bounded. Indeed, from the Holder inequalities < 
^3-e^fl-2 ^^j^ r^^^^ 2004) and (^-3-9^9-2^^ ^ (^'2)?"^(^'3)t^^ one can deduce from 
and 1(213 that 

< (vi/2)r^*3)r' < (2.8) 

This result implies that for G (2,3), (^'e)t is bounded, provided that both {F_i)t and 
{F_2)t are bounded. This condition is assured if s > and Fq is bounded, a condition 
normally required of the forcing, because i^_2 £ F^i/s < Fq/s'^. One may even consider 
a class of forcing for which Fq = oo and F^2 < F^i/s < oo. 

Upper bounds of the above type on dynamical quantities are rather trivial for bounded 
turbulence. However, they are important in the unbounded case, for two reasons. First, 
the scale-selective viscous dissipation allows for the possibility of unbounded growth 
of certain quadratic quantities toward the low wavenumbers. Hence, rigorous bounds on 
dynamical quantities are not as abundant as in the bounded case. Second, analytic studies 
of the nonlinear triple-product transfer function are difficult in unbounded domains. 
In the absence of pointwise estimates for the spectrum, these bounds are particularly 
useful for qualitative estimates of the large-scale distribution of energy. For example, 
pTran (2004)| uses inequality l|2.6|l to argue that the energy spectrum ^'2(fc) should be 
shallower than fc~^, as fc — > 0. 



3. Large-scale energy spectrum 

In this section, it is shown that the physical laws of SQG dynamics admit only large- 
scale energy spectra shallower than k~^. This result is due in part to the fact that 
the simultaneous conservation of '^i and '^2 allows virtually no kinetic energy to get 
transferred toward the low wavenumbers, so that only large-scale kinetic energy spectra 
shallower than A;~^ are possible. 

3.1. Shell-averaged energy spectrum 

For a given wavenumber r, let us denote by S* = S{r) the wavenumber shell between 
k = r/2 and k = 3r/2, i.e. S{r) — {k : r/2 < k < 3r/2}. The shell-averaged energy 
spectrum ^'2(?') over S{r) is defined by 

*2(r) = - / ^2{k)dk. (3.1) 

r J r/2 

In the present case of a doubly periodic domain of size L, the Fourier representation 
of the stream function is ?A(a;) = X]fc ^^Pi*^ ' ^IV'C^)) where k — 2TTL^^{kx,ky) with 
kx and ky being integers not simultaneously zero. Let i/^i^) denote the component of -0 
spectrally supported by S, i.e. i^{S) = J2keS^^P{'^^ ' x}^{k). One has 

sup 1^(5)1 < ^ fc|^^(fc)| <(y.^11 k'lmA < cLr^l^\S), (3.2) 
fees \kes fees / 

where the Cauchy-Schwarz inequality is used, the sum J2kes ^ ~ (cLr)"^ is the number of 
wavevectors in S, c is an absolute constant of order unity and ^'2(6') is the contribution 
to the kinetic energy from S. 
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3.2. Upper hounds for the energy spectrum 

A simple upper bound for \E'2(fc), which is applicable to both the unbounded and bounded 
cases, can be derived from (|2.6ll . In fact, it follows from H2.6II and H3.1|l that 

_ 1 /■3fe/2 

(*2(fc))t = 7/ {^2{>^))tAK<^i-^F_2)tk-\ (3.3) 

« Jk/2 

This bound is supposed to apply to k in the inverse-transfer region. For k in the direct- 
transfer region, lfT7|l yields 

1 ^3fc/2 2 /■3fc/2 

(*2(fc))t = 7/ (*2(«))tdK< - / (vI/3(At))tdK<2Ai-2(i^_^)^^-2, (3 4) 

Jfc/2 Jfe/2 

The upper bound l|3.3|) suggests that dimensional analysis arguments, which predict a 
large-scale k~^ energy spectrum, are not well justified. If a persistent inverse cascade of VPi 
exists (d^'i/dt > 0), then the energy ^'2 ought to acquire a value such that ^'2 < 
In the unbounded case, the large-scale energy spectrum then needs to be strictly shallower 
than k~^, to ensure that the dissipation of does not increase without bound as the 
inverse cascade proceeds toward k ^ 0. On the other hand, if no inverse cascade of \I'i 
exists, then a k~^ energy spectrum with limited extent is possible. If viscous dissipation 
mechanisms with degrees higher than that of the natural dissipation are considered, then 
the upper bounds derived above are not valid. Nevertheless, diminishing energy transfer 
towards the lowest wavenumbers appears to be consistent only with spectra shallower 
than k~^ (for low-wavenumber convergence of the energy integral). The numerical results 
reported in § 4 are well suited to this expectation. 

An upper bound for the large-scale energy spectrum, based on the nonlinear transfer 
term, can be derived for the bounded case. This analysis employs elementary but rigorous 
estimates of the triple-product term. For 3k/2 < s, the evolution of 'i>i{S{k)) is governed 

by 

^<fi{S) = -mS)J{^, (-A)i/V)) - 2m*2(5) 
at 

= ((-A)1/VJ(^,^(^)))-2m*2(5) 

< (|(-A)1/V||VV'||VV'(5)|)- 2^1*2(5) 

< sup |VVX5)|(|(-A)1/2^||V7M) - 2^*2(5) 

X 

< 2cLk^l^'^{S)^2 - 2At*2(5') 

= c^ii~^L^kHl - Aifc*2(fc), (3.5) 

where the second equality is a consequence of ()1.3|l and the second last and last in- 
equalities follow from H3.2|l and the geometric-arithmetic mean inequality, respectively. 
It follows that 

(^2(fc))t < c^^l-^L^k{^l)t. (3.6) 

A notable feature of (|3.t)|) is its dependence on the fluid domain size. The presence of L in 
this upper bound is natural: the upper bound sup^ |VV'(S')|, which is associated with the 
fluid velocity at scales « fc^^, is inherently domain-size dependent. There are no known 
analytic estimates that allow one to derive an upper bound on the nonlinear transfer 
function {i]}{S)J(Tp, (— A)^/^?/;)) in terms of 'intensive quantities' only. This difficulty 
arises not only in the present estimate but also in other analytic estimates of the transfer 
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function. In other words, the nonhnear triple-product term is intrinsically domain-size 
dependent. This problem considerably limits our ability to assess the nonlinear transfer in 
unbounded systems. Finally, it is worth mentioning that although the upper bound H3.6|l 
has a linear dependence on fc, it may be more excessive than the bound {F^2)tk~^ 
derived earlier (even for very low wavenumbers). The reason is that L^k > k^^ and the 
prefactor c^(^2)t '^^y ^'^^ be as optimal as (F-2)t- 



4. Numerical results 

This section reports results from numerical simulations that illustrate the realization of 
large-scale spectra shallower than k""^ . Equation H1.2|) is simulated in a doubly periodic 
square of side 27r, where the forcing /(fc) is nonzero only for those wavevectors k having 
magnitudes lying in the interval K = [59, 61]: 

Here e = 1 is the constant energy injection rate and N is the number of distinct wavenum- 
bers in K. The (constant) injection rate of ^'i is e/s « 1/60, where 1/s k. 1/60 is the 
mean of fc~ over K. This type of forcing was used by [Shepherd (1987)1 |Tran (2004)| 



and Tran & Bowman (2004) in numerical simulations of a large-scale zonal jet on the 
so-called beta-plane and of Navier-Stokes turbulence. The attractive aspect of (I4.1|l . as 
noted in Shepherd (1987) is that it is steady. Dealiased 683^ and 1365^ pseudospectral 



simulations (1024^^ and 2048^ total modes) were performed. Three dissipative forms were 
considered: 2.5 x lO-^Ai/^, -4 x 10"'*(-A)3/2'(/,, and -6 x lO-^A^V + for several 
values of \x. The first case represents the natural dissipation of the SQG dynamics due 
to Ekman pumping, as mentioned earlier. The second case represents thermal (molecu- 
lar) diffusion since (— A)"'^/^-;/) is equivalent to the fluid temperature. The third case — the 
mixed hyperviscous/Ekman dissipation form — is considered in order to demonstrate that 
even slight amounts of Ekman damping will inhibit the formation of an inverse cascade. 



Unlike [Smith tt al. (2002) the case of mechanical friction [oc {—A)^^^^] was not consid- 



ered. The higher resolution was used for the first (natural dissipation) case and the lower 
resolution was used for the second and third cases. All simulations were initialized with 
the spectrum *2(fc) = lO-^Trk/{60'^ + k"^). 

Figure ^ shows the time-averaged steady-state kinetic energy spectrum for the case 
of the natural dissipation term 2.5 x lO~^A-0. The dissipation agents of 5*1 and 5*2 
are, respectively, \['2 (energy) and ^'3. The value of the energy, 0.3333, implies that the 
dissipation of ^'i, averaged in the same period, is 0.01666. This amounts to virtually all 
of the injection rate 1/60. Hence, there exists no inverse cascade of ^1 to the large scales 
and both "ifi and ^"2 are steady. The small-scale energy spectrum scales as fc~^-^, so that 
the spectrum 4*3 (fc) scales as This scaling means that the energy dissipation occurs 

mainly around the forcing region and is consistent with the bound H3.4() . 

Unlike Navier-Stokes turbulence, for which the inverse energy cascade is robust and 
can be simulated at relatively low resolution, it was noticed that no choice for the value 
of ^ at the present resolution could be used to simulate an inverse cascade of ^I^i. It is 
not known whether an inverse cascade of \E'i is realizable at higher resolutions, using 
a smaller value of /i. Nevertheless, this observation suggests that ^'i is 'reluctant' to 
cascade to the large scales, as compared with the more robust inverse energy cascade in 
Navier-Stokes turbulence. 

Figure El shows the kinetic energy spectrum averaged between t = 37.3 and t = 38.7, 
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Figure 1. The time-averaged steady-state energy spectrum ^2(fe) vs. k for the dissipation 

term 2.5 x lO'^Aip. 

for a lower viscous degree. The dissipation agents of "ifi and ^2 are, respectively, ^'3 and 
^4 (enstrophy). The value of V&s is 20, implying that the dissipation of ^'i is 1.6 x 10^^. 
This amounts to about 96% of the injection rate 1/60. The inverse cascade then carries 
only a few percent of the injection of to the large scales. 

The small-scale energy spectrum scales as k~'^-^, so that the enstrophy spectrum '^i{k) 
scales as fc~^-^. Most of the energy dissipation occurs around the forcing region, consistent 
with a 'weak' inverse cascade (one that does not carry virtually all of the injection of 
^1 toward fc = 0; cf. Tran and Bowman 2004, Tran 2004). No direct cascade is possible 
for bounded turbulence in equilibrium or for unbounded turbulence in the presence of a 
weak inverse cascade. 

Similarly, Figure |21 shows the kinetic energy spectrum averaged between t = 15.7 
and t = 16.5 for the mixed dissipation —6 x 10~^A'^ip + fiAip, using three different 
values of fi. When ^ = 0, the dissipation agents of "ifi and are, respectively, ^1/4 
(enstrophy) and ^'5. The value of the enstrophy, 1208, implies that the dissipation of 
\E'i is 1.45 X 10~^, amounting to about 87% of the injection rate 1/60. The small-scale 
energy spectrum scales as k~^, so that ^5{k) scales as k~^. Again, this scahng means 
that most of the energy dissipation occurs around the forcing region and that the inverse 
cascade is weak. We note that as ^ is increased, the inverse cascade becomes increasingly 
weak. We emphasize this behaviour by plotting in Fig. 0] the inverse cascade strength 
r = 1 — 2s{^'^2 + 6 X 10^^^4)/e for six different values of /i. 

Unlike Navier-Stokcs turbulence, for which the enstrophy acquires its near-equilibrium 
value once a discernible inverse-transfer range has formed, the energy in SQG turbulence 
can remain significantly less than its equilibrium value until a very wide inverse-transfer 
range has developed. For example, for a one-decade Navier-Stokes inverse-transfer range 
(achievable in numerical simulations), the enstrophy acquires 95% of its projected equi- 
librium value (calculated with a energy spectrum extrapolated to fc = 0). On the 
other hand, for a one-decade SQG inverse-transfer range, the energy acquires only 66% of 
its projected equilibrium value (calculated with a fc~°'^ energy spectrum extrapolated to 



C. V. Tran and J. C. Bowman 

10- 




10- 



10" 



10' 



10= 



103 



Figure 2. The quasisteady energy spectrum '^2{k) vs. k averaged between t — 37.3 and 
t = 38.7 for the dissipation term -4 x 10~^{-A)^/^ip. 
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Figure 3. The quasisteady energy spectrum ^'2(fc) vs. k averaged between t = 15.67 and 
t = 16.52 for the dissipation term —6 x 10~^A'^i/> + iiAip, using three different values of fi. 



k — 0, as realized in the present simulations; cf. the ii — case of Figure |3Jl. This means 
that one needs a considerably wider inverse-transfer region for SQG turbulence than for 
Navier-Stokes turbulence, in order to approach a quasi-steady state. This problem is in 
addition to the resolution limitations at the small scales for both cases. 

Due to the steep spectrum in the inverse-transfer region, the energy in the /i = case 
of Figure |31 has not acquired a value considerably close to its equilibrium value. This 
means that the system is still well within the transient phase, However, the dissipation 
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Figure 4. The decay of the inverse cascade strength r for the dissipation term 

as IX is increased. 

of ^'i (proportional to the enstrophy) cannot grow considerably (without significant 
change to the existing spectrum), because of the high degree of viscosity, which makes 
the dissipation of relatively insensitive to growth of the large-scale energy. 

5. Conclusion and discussion 

In this paper, the kinetic energy density of SQG turbulence and its large-scale spectrum 
have been studied. For the unbounded case, upper bounds arc derived for the time means 
of the kinetic energy density and of the large-scale energy spectrum, averaged over a 
narrow window of wavenumbers. Another result is an upper bound on the the time mean 
of the large-scale energy spcctriim, which is derived for the bounded case. Numerical 
results confirming the predicted slopes of the large-scale energy spectrum are presented 
and discussed. 

An important feature in SQG turbulence that gives rise to the rigorous upper bound 
on the time mean of the kinetic energy density in the unbounded case is that the kinetic 
energy is the dissipation agent of the inverse-cascading candidate '^i. This fact is due to 
the hypoviscous nature of the dissipation operator (— A)"'^/^, a natural physical dissipation 
mechanism of SQG dynamics (Ohkitani 1997; Constantin 2002; Tran 2004). If (-A)i/2 
is replaced by an operator of the form (—A)'', where rj > 1/2, then the simple analysis 
of Section 2 fails to show that the time mean of the energy density {^2)t is bounded, 
although it may remain so for low degrees of viscosity rj. The reason is that the amount of 
energy getting transferred to wavenumbers lower than a given wavenumber k decreases 
at least as rapidly as /c, so that the spectral dissipation rate oc fc^'', a consequence of the 
dissipation operator (— A)"*, may be sufficiently strong to balance the diminishing inverse 
energy transfer and keep the energy from growing unbounded. 

Numerical simulations of SQG turbulence were performed, using the natural dissipation 
operator (— A)^/^ and two viscous operators A and (— A)"^/^. The results show large-scale 
energy spectra shallower than k~^, consistent with the theoretical prediction. 
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There have been attempts to explain, within the context of SQG turbulence (Con- 
stantin 2002; Tung & Orlando 2003), the kinetic energy spectra observed in the laboratory 



experiment of Baroud (2002) and in the atmosphere. In the former case, the turbulence in 
a rotating tank is driven at a sufficiently small scale to allow for a wide inverse-transferring 
range. A spectrum extending over nearly two wavenumber decades lower than the 
forcing wavenumber is observed. In the latter case, a k~^/^ spectrum is observed in the 
mesoscales (see Frisch 1995 and Tung & Orlando 2003 and references therein), which 
correspond to wavenumbers higher (lower) than the forcing wavenumber if the energy 
released from baroclinic instability (thunderstorms) is considered to be the driving force. 



The —2 power-law scaling observed in [Baroud (2002) for the wavenumber range lower 



than the forcing wavenumber is excessively steeper than the permissible scalings derived 
in this work. The —5/3 slope in the atmosphere is either steeper (if considered to be on 
the wavenumber range lower than the forcing wavenumber) or shallower (if considered 
to be on the wavenumber range higher than the forcing wavenumber) than the permis- 
sible slopes. According to the present analysis, these data cannot be attributed to SQG 
turbulence. 

We would like to thank two anonymous referees for their comments, which were help- 
ful in improving this manuscript. This work was funded by a Pacific Institute for the 
Mathematical Sciences Postdoctoral Fellowship, an Alexander von Humboldt Research 
Fellowship, and the Natural Sciences and Engineering Research Council of Canada. 
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